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In the framework of risk management, for the study of the sensitivity of pricing and 
hedging in stochastic financial models to changes of parameters and to perturbations 
of the stock prices, we propose an error calculus which is an extension of the Malliavin 
calculus based on Dirichlet forms. Although useful also in physics, this error calculus is 
well adapted to stochastic analysis and seems to be the best practicable in finance. 

This tool is explained here intuitively and with some simple examples. 
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1. INTRODUCTION 

Once a model is chosen to price contingent claims and to hedge a position, the 
actual questions obviously are : what is the exposure to errors on the model and to 
changes in the market ? This risk assessment is usually done in terms of sensitivity 
of the portfolio to variations of numerical financial quantities and parameters of the 
model. Now the theory of Dirichlet forms allows to extend this sensitivity calculus 
to perturbations of functional quantities like stochastic processes. The method 
is an extension of Malliavin calculus. The errors are thought to be infinitesimal 
random quantities with biases, variances and covariances. 

Among the promising clues of application of error calculus to finance, let us 
mention some directions : to manage the precision of numerical methods used to 
implement the stochastic theory ; to obtain new integration by parts formulas to 
speed up Monte Carlo simulations ; to study how depends the solution of an ODE 
or of an SDE on a functional coefficient like the level-type volatility. We sketch 
these two last directions at the end of the article. 

Our presentation starts from the basic ideas and goes until examples of com- 
pletely tractable computations. 

2. PRESENTATION OF THE METHOD 



Several technics are available to represent errors mathematically and to compute 
them with formulas. The method of Dirichlet forms is based on some major his- 
torical ideas which are the easiest way to penetrate it. 



2.1. Propagation of errors : error calculus a la Gauss 

After his argument showing the importance of the normal law ( Theoria motus cor- 
porum coelestium 1809), Gauss was interested in the propagation of errors (Theoria 
Combinationis 1821). Given a quantity U = F(Vi, V 2 , ■ ■ ■) function of other erro- 
neous quantities Vi, V 2 , . . . he states the problem of computing the quadratic error 
to fear on U knowing the quadratic errors of , of , . . . on V 1: V 2 , ■ ■ . , these errors 
being supposed small and independent. His answer is the following formula 




he gives also the covariance between the error on F and the error of an other 
function of the V$'s. 

Formula (2.1) possesses a property which makes it highly better, in several 
questions, than other formulas used here and there in textbooks during the 19th 
and 20th centuries. It is a coherence property. By lack of place we refer to Bouleau 
(2001) §1 for the reason of this coherence property. 

In the calculus a la Gauss the errors on Vi, V 2 , . . . are not necessarily supposed 
to be independent nor constant, they can depend on V±, V 2 , . . . : Let be given a field 
of symmetric positive matrices (o^-^i, v 2 , . . .)) on IR d representing the conditional 
variances and covariances of the errors on Vi, V 2 , . . . given the values v 1, v 2, . . . of 
Vi, V 2 , . . . then the error on U — F(Vi, V 2 , . . .) is 

(2.2) a\ = |f v 2 , . . -)§^-(vi, v 2 , . . .)a ij (v 1 , v 2 , . . .) 

which depends solely on F as mapping. 

2.2. First order and second order calculus 

The following remark, although very simple, is important to understand the 
role of the error calculus a la Gauss that will be used in the sequel in the extended 
form allowed by Dirichlet forms. 

Let us start with a quantity x with a small centred error eY, on which acts a 
non-linear regular function /. Thus we have at the beginning a random variable 
written x + eY, it has no bias (centred at the true value x) and its variance is 
e 2 a 2 Y . 

After having applied the function /, using Taylor formula shows that the error 
is no more centred. The bias has the same order of magnitude as the variance. 
Then applying new regular non-linear functions f n gives a transport formula which 
shows how errors propagate : biases and variances keep permanently the same 
order of magnitude 

bias n+ i = bias n /; +1 (x„) + ivariance„/;' +1 (a; n ) + e 3 0(l) 
variance n+ i = variance n /^ + i(a;„) + e 3 0(l) 

(it could be easily extended to applications from 1R P to JR g , for the general formulas 
on the bias and the variance of the error under regular mappings see Bouleau and 
Hirsch (1991) chapter I §6 corollaries 6.1.3 and 6.1.4). 



We see that the calculus on the biases is a second order calculus involving the 
variance. Instead, the calculus on the variances is a first order calculus not involv- 
ing the biases. Surprisingly, the calculus on the second order moments of errors is 
indeed simpler than the calculus on the first moments. Thus, the error calculus 
on the variances appears to be necessarily the first step in an analysis of errors 
propagation based on differential methods and supposing small errors. 

2.3. Extended error calculus using Dirichlet forms 

The error calculus of Gauss has the limitation that it has no mean of extension. 
If the error on (Vf, V 2 , V3) is known it gives the error on any differentiable function 
of (V 1 ,V 2 ,V 3 ) but that's all. 

Now, in the usual probabilistic situations where a sequence of quantities X ± , X 2 , . . 
is given and where the errors are known on the regular functions of a finite number 
of them, we would like to deduce the error on a function of an infinite number of 
the Xj's or at least on some such functions. 

It is actually possible to reinforce this error calculus giving it a powerful ex- 
tension tool and preserving the coherence property. In addition, it will give us the 
comfortable feature to handle Lipschitz functions as well. 

For this we come back to the idea that the erroneous quantities are themselves 
random, as Gauss had supposed for his proof of the 'law of errors', say defined on 
(fl,A, IP). The quadratic error on a random variable X is then itself a random 
variable that we will denote by V\X\. Intuitively we still suppose the errors are 
infinitely small although this doesn't appear in the notation. It is as if we had an 
infinitely small unit to measure errors fixed in the whole problem. The extension 
tool is the following, we assume that if X n — > X in L 2 (Q,A, P) and if the error 
T[X m — X n ] on X rn — X n can be made as small as we want in A, IP) for m, n 

large enough, then the error T[X n — X] on X n — X goes to zero in L 1 . 

This can be axiomatized as follows : we call error structure a probability space 
equipped with a local Dirichlet form possessing a carre du champ. 
Definition 2.1. An error structure is a term (Jl, A, P, ID, T) where (Q,A, P) is 
a probability space, satisfying the four properties : 

1. ) E) is a dense subvector space of L 2 (fl, A,JP). 

2. ) r is a positive symmetric bilinear map from DxD into £ X (P) fulfilling the 
functional calculus of class C 1 n Lip, what means that if u E JD m and v G P ra for 
F and G of class C 1 and Lipschitz from P m [resp. JR n ] into JR, one has Fou G P 
and GodgD and 



3.) the bilinear form £[f, g] = JET[f,g] is closed, i.e. P is complete under the 



With this definition, the form S defined at point 3.) is a Dirichlet form. This 
notion has been introduced by A. Beurling and J. Deny as a tool in potential theory, 
see Beurling and Deny (1958-59), and also Fukushima, Oshima and Takeda (1994). 



T[Fou,Gov] = J2 t , J F;(u)G' J (v)T[u l ,v J ] JP-p.s.. 




4.) 1 G P and V[\, 1] = 0. 
We always write S[f] for S[f, f] and T[f] for r[/, /]. 



The operator T is the carre du champ or squared field operator associated with £, 
it has been studied by several authors in more general contexts, see Dellacherie and 
Meyer (1987), Bouleau and Hirsch (1991). Here we refer to T as the quadratic error 
operator of the error structure. Its intuitive meaning is the conditional variance of 
the error. 

Example 2.1. A simple example of error structure is the term 

(IR.SQR),™,// 1 ^)^) 

where m is the normal law N(0, 1) and 

if 1 (to) — {/ G L 2 (m) : /' in the sense of distributions G L 2 (m)} 

with 7 [/] = f' 2 for / G if 1 (to). This structure is associated to the real valued 
Ornstein-Uhlenbeck process. It models an erroneous quantity, say X, with nor- 
mal law whose error does not depends on the value of X. Instead, the operator 
l[f]( x ) = f' 2 {x)x 2 would model an error proportional to X, that is, in the sense 
of physicists, a constant proportional error. The intuitive relation giving the in- 
terpretation of the quadratic error operator 7 is 

7 [/](x) = E[(error on f{X)f\X = x\. 

2.4. How proceeds an error calculation 



Let us suppose we are drawing a triangle with a graduated rule and a protractor: 
we take the polar angle of OA say 9±, we put OA = £1, then we take the angle 
(OA, AB) say 9 2 , and we put AB = £ 2 . 

1) Choose hypotheses on errors 

£±, £2 and 6>i, 9 2 and their errors can be modeled by the following probability 
space and operator : : 

((0, Lf x (0, tt) 2 , £?((0, L) 2 x (0, tt) 2 ), f f ID, V) 

where JD = {/ G L 2 (f f ^) : f , f , J£, J£ G L 2 (f f f^)} and 

It is easily checked that assumptions 1) 2) 3) 4) of definition 2.1 are fulfilled. 

2) Compute the errors on significant quantities using the functional calculus on 

r 

For the coordinates of the point B for example we have : 

X B = £ 1 cos9 1 +£ 2 cos(9 1 + 9 2 ) Y B = l x sin 9 X + £ 2 sin{9 1 + 9 2 ) 

T[X B ] = £\ + £ x £ 2 (cos # 2 + 2 sin 9 1 sin(0i+0 2 )) + £ 2 {l + 2 sin 2 (^ +0 2 )) 

r [Tb] = ^ + £i£ 2 (cos # 2 + 2 cos 9 l cos (^ + 9 2 ) ) + £ 2 ( 1 + 2 cos 2 (9 1 +9 2 )) 

T[X B , Y B ] = -£ x £ 2 sin(29 1 +9 2 ) - £\ sm(29 1 +29 2 ) 

Then, according to the problem, we can for example compute the covariance of 
the errors on the area and on the perimeter of the triangle, etc., or obtain that the 
proportional error on the area T [area] /area 2 is maximal for 9 2 = ir/2, etc. 



Remark 2.1. If we limit our investigation to variances of the errors, that is to 
computation with V, then the choice of the a priori laws is not so crucial as it could 
be thought because these computations are done almost surely (using property 2 
of definition 2.1. If we are, instead, interested also in biases, then the a priori laws 
are precisely relevant. Biases are represented by an operator which is the genera- 
tor of the semi-group canonically associated with the error structure (see Bouleau 
and Hirsch (1991) chapter I). In example 2.1 it is A [«](#) = \u"(x) — \xu'(x). If 
we change the probability measure m into f.m (/ regular), 7 being unchanged, 
then the operator A becomes A[u] — A[u] + 7^7 [f,u]- The operator u — > jjj[f,u] 
is first order. Absolutely continuous changes of the probability measure m corre- 
spond to changes of the drift of the bias operator A, a variant of Girsanov theorem. 

2.5. Comparison of approaches 

Before looking at the infinite dimensional examples needed in finance, let us 
try to give an outlook over the different approaches to error calculus. 

At the extreme right-hand side of the table we have the usual probability 
calculus in which the errors are random variables. The knowledge of the joint laws 
of the quantities and their errors is supposed to be yielded by statistical methods. 
The errors are finite, the propagation of the errors needs computation of image 
probability laws. 



deterministic 
approach 


probabilistic approaches 


Sensitivity 
calculus: 
derivation with 
respect to the 
parameters of 
the model 


Extended error 
Dirichk 
first order calcu- 
lus only dealing 
with variances 


calculus using 
>t forms 
second order cal- 
culus with vari- 
ances and biases 


Probability 
theory 


infinitesimal errors 


finite errors 



Table 2.1: Main classes of error calculi 



At the extreme left-hand side the usual sensitivity calculus consists of comput- 
ing derivatives with respect to parameters. Let us remark that it applies also to 
functional coefficients using Frechet or Gateaux derivatives. 

Between these two purely probabilistic and purely deterministic approaches 
lies the extended error calculus based on Dirichlet forms. It supposes the errors 
infinitely small but takes in account some features of the probabilistic approach 
allowing to put the computations and the arguments inside a powerful mathemati- 
cal theory: the theory of Dirichlet forms. In the same framework can be performed 
either a first order calculus on variances which is simple and significant enough for 
most applications or a second order calculus dealing with both variances and biases. 

2.6. Main features of the method 

As above in the finite dimensional case of the triangle, the construction of an 
error structure on an infinite dimensional stochastic model is done in two steps 



1) If there are, as usually, deterministic parameters which can be erroneous 
or with respect to which a sensitivity is wished, these parameters have to be 
randomized with a priori laws. 

2) Errors operators must be chosen to act on random quantities (initially ran- 
dom or randomized parameters) in order to describe errors, in such a way that we 
obtain mathematically an error structure. 

Several properties of error structures make it easier such a construction. 

1) The operation of taking the image of an error structure by a mapping is quite 
natural and gives an error structure as soon as the mapping, even non injective, 
satisfies some rather weak conditions. In particular if (ft, A, P, ID, T) is an error 
structure and if X is a random variable with values in 1R d whose components are 
in E>, (]R d ,i3(IR d ),Px,lDx,rx) is an error structure where JP X is the law of X 
and 

E> x = {feL 2 (F x ):foXED} 

r x [f} = m\r\foX]\x=x], feJD. 

2) If / G D and F is Lipschitz from 1R to IR then Fof g ID and T[Fof] < kT[f]. 
For example the structure of example 1 (H, £>(H), m, i? 1 (m), 7) possesses an image 
by the map x — > | sin y / T-T~|xj| which is an error structure on [0, 1]. Such a use 
of non injective functions is tricky in the deterministic sensitivity calculus. More 
generally if F is a contraction from JR d into IR in the following sense 

\F{x)-F{y)\<Yf i=l \x i -y i \ 

then for f u / 2 , . . . , f d G ID one has F(f 1} f 2 , . . . , f d ) G ID and 

r[F(/ 1 ,/ 2 ,...,/ d )]UEi 1 r[/ i ]i 

This property allows to consider more general images with values in metric spaces 
as soon as a suitable density property is preserved, see Bouleau-Hirsch (1991) 
chapter V §1.3 p 197. 

3) The product of two or countably many error structures is an error structure. 
It is the mathematical expression of the independence of the random variables and 
the non-correlation of the errors. By this way error structures on infinite dimen- 
sional spaces are easily obtained, e.g. on the Wiener space, as we will see in the 
next part, or on the general Poisson space or other spaces of stochastic processes, 
see Bouleau and Hirsch (1991), Ma and Roeckner (1992), Bouleau (1995). 

For later reference we give the following statement. 
THEOREM 2.1. Product structures 

Let S n = (Q n , J^n, m n , D n , T n ), n > be error structures. 

The term S = (fi, J 7 , m, ID, T) defined below is an error structure denoted S = 
Il^Li S n and called the product structure of the S n : 

(n, F, m) = (TXLo ^n, ®^=o I\n=o m n ) 

ID = {/ G L?{m) : Vn, for m-a.e. uj = (u , ui, . . .) 

the function x — > f(u , . . . , u; n _i, x, u n+ i, . . .) G ID n 
and /E„r n [/] dm < +00} 

and for fe ID F[f] = E„ T n [/]. 



Thanks to these properties, is possible the construction of a variety of error 
structures on a given probabilistic model. Now for a rational treatment of a practi- 
cal case these error hypotheses should be obtained by statistical methods. This is 
connected with the Fisher information theory, see Bouleau (2001). Anyhow, these 
statistical methods are not yet sufficiently studied to be exposed here, especially 
in the infinite dimensional case we have to use in finance. Thus we limit ourselves 
to error computations with a priori errors chosen the most likely we can. We will 
see that it is significant already. 

3. ERROR STRUCTURES ON THE WIENER SPACE 

Let us first recall the classical construction of the Brownian motion using the 
Wiener integral. 

3.1. The Wiener space as Gaussian product space 

Since we aim here at applications we will consider only the case where a mea- 
sured space (E, S, //) is given which is either (IR+, B(JR + ), dt) or ([0, 1], £>([0, 1]), dt) 
and a one-dimensional Brownian motion (for the abstract Wiener space setting see 
Bouleau and Hirsch (1991)). 

Let (xn) be an orthonormal basis of L 2 (E,£,fi) and let (g n ) be a sequence of 
i.i.d. reduced Gaussian variables defined on a probability space (f2, A, IP). To each 
/ G L 2 (E,£,/j,) we associate /(/) G L 2 (tt,A,1P) by 

Hf) =En < f,Xn> 9n- 

then / is an isometric homomorphism from the Hilbert space L 2 (E,£,fi) into the 
Hilbert space L 2 (fl, A, IP). If / and g are orthogonal in L 2 (E, £, /i), 1(f) and 1(g) 
are independent Gaussian random variables and putting 

(3.1) B t = En < l[o,t], Xn >9n (t G [0, 1] or t G R+) 

defines a Gaussian stochastic process which is easily shown to be a standard Brow- 
nian motion. By extending the case where / is a step function, the random variable 
/(/) is denoted by / f(s) dB s and defines the Wiener integral of /. In this con- 
struction we can suppose the space (Q, A, P) be a product space: 

(£l,A,X>) = (P,£(P),m) N m = N(0,l) 

and the g n 's be the coordinate mappings. Thus uj = (u> , . . . , u> n , . . .) and g n (uS) = 

UJ n . 

By the functional calculus, as soon as ID and V define an error structure on 
(Q, A, P), say (Q, A, P, E), T) for which the g n 's are in E), this structure is deter- 
mined by the quantities 

(3.2) r[J7(s) dB s ] f G A) dense in L 2 [0, 1] (resp. L 2 [P+]) 

because it follows that if F G C l n Lip(R k ) 

T[F(ff 1 dB,...,ff k dB)] = 

Et=i m fidB, . . .) F>(f f ± dB, . . .)T[f f,dU. f fjdB] 



and the random variables F(J fidB, . . . , / f k dB) for F G C 1 nLip(WL k ) and /; G £> 
are a dense subspace of L 2 (P). 



3.2. The Ornstein-Uhlenbeck structure 



Taking T[J f(s) dB s ] = \\f\\ 2 L 2 and D = L 2 [0, 1] gives a closable structure 
which is of the form 

(Q, A, P, ID, r) = Un= m, B(B),m, d n , 7 „) 

where each factor (P, £>(P), m, d„, 7n ) is here a copy of (P, £>(P), m, if 1 (m), 7) 
the Ornstein-Uhlenbeck structure of example 2.1. This error structure is induced 
by the following perturbation of the Brownian path : 

(3.3) u(t) — ► V(F°u(t) + \jl-e- 6 u(t) 

where Co is an independent Brownian motion and 6 a vanishing parameter. 



3.3. Finer structures of product type 



Taking T[J f(s) dB s ] = \\f'\\ 2 L2 and D = C^O, 1] or more generally for q G IN* 
(or even q G P*J T[J f(s) dB s ] = \\f {q) \\ 2 L 2 and D = C[0, 1] give also closable 
structures. For a suitable choice of the basis (% n ) they are also of the form 

(Q, A, P, ID, r) = Un=o((^ B{JR),m, d n , 7n ). 

3.4. Error structures of generalised Mehler type 

Let p t be a strongly continuous symmetric contraction semi-group on L 2 (P + , dt) 
with generator (a, T>a) (a being a non necessarily local operator), let us consider the 
associated closed positive quadratic form (e, V(\f^a)) defined by e[f] = \\ \J — a/|| 2 
(a non necessarily Dirichlet form), then the structure on the Wiener space induced 
by the formula 

r[f °°f(s)dB s ]=e[f], feV(^a)) 

is closable and thus defines an error structure. 

It corresponds to the semigroup P t on L 2 (f2, A, P) given by 



P f F = E 



(/•CO roc J _ 

J o (pt/2l[o,.])(u)dB u + J o {^I- Pt l [QA ){v)dB v 



where \/l — p t is the positive square root of the positive operator I — p t on 
L 2 (]R + ,dt) and B is an auxiliary independent Brownian motion. 
Remark 3.1. When a financial model is studied by means of a development in 
series with respect to a small random change in the coefficients (like small noise 
expansion of a stochastic volatility) it is possible to induce from the perturbation 
an error structure which manages the variances and the biases at the limit when 
the perturbation is infinitely small. We cannot describe the details here of this 
standard method. That yields often error structures outside the class of general- 
ized Mehler type, but still defined by the quantities (3.2). 



3.5. The gradient operator and the derivative 



In any error structure whose space ID is separable, we can define a gradient 
operator D on ID with values in L 2 (P, H) where H is an auxiliary Hilbert space: 
D is a continuous application from ID into L 2 (1P, H) such that 

1) W, V G ID < DU, DV > H = T[U, V] 

2) VF G C 1 n Lip(WL d ), VIgD D(F o X) = £ti F/ o XM* P-a.s. 
Then if £/, V G ID n L°° which is an algebra, it holds D(UV) = DU.V + U.DV. 
For example in the case of the Ornstein-Uhlenbeck structure, taking H = L 2 (IR + ) 
gives 

. V/i G L 2 (JR+) D ou [J h{s) dB s ] = h 

. with suitable hypotheses on the adapted processe H t (see Nualart (1995)) 

D ou [JH s dB s }(t) =H t + !{D ou H s ){t) dB s 
.HUE ID ou the Clark formula 

U = WJ + / °° m[D m U(t)\F t ] dB t . 
Now a slight variant of the gradient operator, the notion of 'derivative* , is useful 
when computing errors on solutions of stochastic differential equations thanks to 
the tool of Ito's formula (this notion has been used and studied by Feyel and la 
Pradelle (1989)). 

Definition 3.1. Let (B t )t>o be an auxiliary independent Brownian motion. For 
U G ID the derivative is a random variable depending on u and Cj defined by 

U* = JS°(D m U){u,t) dB t . 

From the properties of the gradient one gets 
.T[U}= IE[f/# 2 ] 

. Fox FeC 1 D Lip (F o U)* = F' o U . U* 
(J H s dB s )* = JH s dB s + J Hf dB s . 

3.6. The weighted Ornstein-Uhlenbeck case 

Its meaning for financial models is to consider non necessarly time translation 
invariant perturbations of the underlying stock price. It is a special case of the 
generalised Mehler type: 

T[J f(t)dB t ] = J a(t)f 2 (t)dt f G V(R+) a measurable > 

. V/i G L 2 (IR + , (1 + a)dt) D[j h(s) dB s ] = y/a(jjh 

. (J h(s) dB s )# = J sj^t) h(t) dB t 
. with suitable hypotheses on the adapted processe H t 
D[JH s dB s }(t) = ^At)H t + !{DH s ){t) dB s 

(J H s dB s )# = J H s ^{s)dB s + / Hf dB s 
. If U G ID n ID 0U DU = ^D ou U 
The generator (A, T>A) of this structure can easily be seen to verify 

A[J fdB} = -J a(s)f(s)dB s f G V(JR+) 

which permits (see formula (JbM)) in the concluding remarks) to compute A on a 
dense part of VA. 



In the sequel, we focuse on the Ornstein-Uhlenbeck case, but the two following 
lemmas are also valid in the weighted Ornstein-Uhlenbeck case, and part III and 
IV extend to that case with only minor changes. 

Lemma 3.1. The conditional expectation operators JE[. |jF t ] are orthogonal projec- 
tors in ID on errors sub-structures (closed sub-vector-spaces of J) stable by Lipschitz 
functions). 

Lemma 3.2. Under the same hypotheses, let T t be defined from T by 

T t [(Jf(s) dB s )]=T[(Jl m f(s) dB s )} 
and let U — > U #t the derivation operator associated with T t , then for U G ID: 

(m[u\r t ])# = m[u**\r t ]. 

4. APPLICATION TO FINANCIAL MODELS 



4.1. The Black-Scholes case 



Notation 

The interest rate for the bond is constant, the asset (S t )t>o is modeled as the 
solution of the equation dSt = St(fidt + adB t ). For a European option of the form 
/(S't), T fixed deterministic time (see Lamberton and Lapeyre (1997)), the value 
at time t G [0, T] of the option is V t — F(t, S t , a, r) with 

-r(T-t) f f , < r -^)(T-t)+ayVT-t\ e 2 



(4.1) F(t, x, a, r) = '*> / /(xe^ W -*J wr- t) - dy 

Jn V27T 



If / is Borel with linear growth, the function F is C 1 in t G [0, T[, C 2 and Lipschitz 
in x g]0, oof, let us put 

delta* = |f (t, S t , a, r) gamma t = 0(t, S t , a, r) 

F satisfies the equation ^ + ^-^r + r ^ff — rF = 0. 
Hypotheses 

Our choice is governed by an aim of simplicity. 

a) The error on (B t ) t > is represented by the Ornstein-Uhlenbeck error struc- 
ture. 

b) The errors on the initial value So, on the volatility a, on the rate r are 
'constant proportional errors' in the sense of physicists : 

iWo)]= ns )si 

rK(r)]= e' 2 (r)r 2 

c) We chose a priori laws : lognormal laws on S and a, an exponential law on 

r. 

d) We suppose (B t ) t > and the randomized quantities are independent and 
their errors uncorrelated. (In a more complete study, these independence and 



uncorrelation assumptions would have to be relaxed, in particular to express links 
between errors on the asset (St) and on the volatility a). 

In other words, the error on a regular function F((B t ) t >o, So, a, r) will be rep- 
resented by the product error structure i.e. 

T[F((B t )t>o, So, a, r)) = Y ou [F(., S , a, r)\ + F' 2 o S 2 + F l2 a 2 + F' r 2 r 2 

where T ou is the Ornstein-Uhlenbeck quadratic error operator. 

Actually, the theory tells us that hedging and pricing formulas do not involve 
the drift coefficient So we may take /x = r, i.e. we work under the probability 
IP such that S t = e~ rt S t , the discounted stock price, is a martingale. Since S t = 

Soe aBt+ ^-^ )l we have 

V[S t \ = S 2 {a 2 J* a(s)ds + (B t - otfa 2 + t 2 }. 
4.2. Errors on the value and the hedge of a European option 

Let us consider an option of the form /(St) where / is Lipschitz. 
By the independence hypothesis, the errors on B, So, o~, r can be managed 
separately. Let us denote T B , r , IV, T r the corresponding quadratic operators. 

a) Error on the value of the option 

The value of the option is V t = F(t, S t , a, r) with F given by (14.11) 
al) Error due to B 

B being present only in S t , we have T B [V t ] = (§f(t, S t ,o-,r)) 2 T B [S t ] so 

U9] r B [v t ]= delta t 2 r B [^] 

{ ' r B [V s ,V t }= delta s delta t T B [S s ,S t ] 

with T B [S s ,S t ] = S s S t a 2 (sAt). 

Proposition 4.1. If f is Lipschitz, V t is in ID B and when t] T 
Vt = F(t, St, o~, r) — > f(Sr) in!D B and P — a.s. 
T B [V t ] = (de\ta t ) 2 T B [S t ] -> f' 2 (S T )T B [S T ] in L 1 and IP - a.s. 

proof. Let us suppose first / G C 1 PI Lip. By the relation 

V t = ne- r(T - t) f(S T )\F t ] 

it follows that V t — > f(Sr) in L p 1 < p < oo and a.s. 

A computation that we shall do in a more general framework later, and that 
we do not repeat here, gives 

V* = e- r ^nf{S T )S T \T t ]oB t 

thus 

V t * -> f(S T )S T aB T in L 2 (JP,L 2 ((l,P)) 
and thanks /(St)* = f (St)St<jBt we obtain 

V t -> f(S T ) in JD B and P - a.s. 



and 

T B [V t ] = e-^-^infiS^ST^faH - f' 2 (S T )T B [S T } 

in L l and IP-a.s. □ 
The case / only Lipschitz comes from a special property of the one-dimentional 

functional calculus in error structures (see Bouleau and Hirsch (1991) chapter III 

prop. 2.1.5), the preceding argument still remains valid. 

a2) Error due to a. 

We suppose here / G C 1 R Lip. As V t = F(t, S t , o~, r) 

TaWt] = {F^((t, S t , a,r)^ + F^(t, S t , a, r)} 2 a 2 
and the computation can be done using the integral representation (|4.1|) . puting 

F(t, x, a) = e r(T -' } F(t, xe"^^, a, r) 

and remarking that by (|4.1|) we have 

dF | <t 2 x 2 9 2 F _ n a£ _ _2(T-t) a£ dF _ / T + \ 2d 2 F 

dt ^ 2 dx 2 ~ V da ~ a dt da ~ V J V OX dx' 2 

and we obtain 

r a [V t ) = {(T — t)(j5 2 gamma t + S t (.B t - (j^deltaj 2 a 2 

(4.3) = {vega 4 + S t (B t - at)delta t } 2 a 2 

One gets immediately, for example, the well-known fact that for two European op- 
tions of payoffs /(i) (St) and /(2)(Sr), an option with payoff ai/(i)(Sr) +02/(2) (St) 
would have a value Vo at t = insensitive to a, i.e. IV [Vb] = 0, as soon as 

(!) 1 ( 2 ) n 

ai gamma + 02 gammajj = (J. 
a3>) Error due to r. 
We have similarly 

T r [V t ] = {F^((t, S t , *,r)-^ + S t , a, r)} 2 r 2 

thus 

(4.4) T r [V t ] = {tStdelt&t + rhoj 2 r 2 . 

As a consequence, given several options of payoffs /(^(St), i — 1, . . . ,k, the option 
of payoff J2i a if(i){ST) has a value at t — insensitive to both a and r (i.e. IV [Vo] = 
r r [Vo] = 0) if the vector a = (<2j) is orthogonal to the two vectors (gamma|, ! ') and 
(rho«). 

b) Error on the hedging portfolio 

Here we limit ourselves to the error due to (B t ). We suppose / and /' in 
C 1 fl Lip. The hedging equation is 

e~ rt F(t, Si, a, r) = F(0, S , a, r) + J* H s dS s 

where the adapted process H t is the quantity of stock in the portfolio : 

BF 1 
H t = delta* = —{t,S t ,a,r) = e- r ^m[f(S T )S T \T t }-. 



By the same method as for V t we obtain 



T B [H t }= (gamma,) 2 ^] 
r B [H s ,H t ]= gamma s gamma,r B [5 s , S t ] 

Proposition 4.2. // /, /' e C 1 n Lip, then H t e ID and ast]T 

H t — > /'(St) in 1D b and a.s. 
T B [H t ] -> /" 2 (S r )r B [S r ] in L^P) and a.s. 



Remark 4.1. These results show that the Greeks introduced by practioners have 
a direct sense as sensitivity of the value V t and of the hedging H t to perturbations. 
This is of course not surprising, the method makes more precise the correlations 
of errors. It gives also a tool to study the absolute continuity of joint laws as we 
explain now. 

The preceding computations show easily that in the Black-Scholes model, if U\ 
and U2 are two random variables taken among the following quantities defined at 
a fixed instant t: S t , Vt(/i), V t (f 2 ), #t(/i), H t (f 2 ), then the matrix Tp^Uj] is 
singular: the errors on these quantities are linked. This comes from the fact that 
the law of e.g. the pair (V t (fi), V t (f2)) is carried by the A-parametrized curve: 

y= exp-r(T-*)iWi(A) 
x= exp-r(T -t)P T _ t f 2 (\) 

where (P t ) is the transition semigroup of (S t ) . The same phenomenon happens in 
any more general Markovian model. 

On the contrary the random quantities involving several different instants have 
generally nondinked errors. Thus for example if U\ = St and £/ 2 = J T e~ s H s S s ds 
(discounted immobilization of the portfolio) the matrix T[Ui,Uj] is a.s. regular as 
soon as / is not constant, hence, by the absolute continuity criterion (Bouleau and 
Hirsch (1986) or Nualart (1995) thm 2.1.2) the law of the pair (S T , Jo e' s H s S s ds) 
possesses a density. 

c) More general errors on (B t ) 

The relations (4.2) (4.3) (4.4) (4.5) still hold in the weighted Ornstein-Uhlen- 
beck case, and also, with suitable hypotheses, if we consider more general error 
structures on the Wiener space. Let us consider, as mentioned above, a structure 
induced by a closed positive quadratic form e on L 2 (IR + , dt) with 

T B [JfdB]=e[f} 

for / in the domain of e with, for example, 

i) e[f] = JS°(f( 8 ) - f(t)) 2 (3(s)(3(t) dsdt 

ii) eif]=fo(f (q) (s)) 2 ds 

where is the fractional derivative of order q, 
then the formulas 

T B [V t ]= (delta*) 2 T B [S t ] 
T B [H t ]= (gamma,) 2 !^] 



remain valid as soon as St G ID i.e. 

in case i) if (3 G L l (M + ,dt) and T B [S t ] = S 2 a 2 2 f t °° (3(s)ds j£ (3(s)ds. 
in case ii) if q G (0, \) and T B [S t ] = S 2 a 2 ££= 1 ^g 1 . 
Now in the case 

iii) e[/]=/o~(E?=iOi(*)/ ( W<k 

we do not have anymore l[o,t] G dom(e), hence B t doesn't belong to ID. Such error 
structures are more convenient to model errors on processes with finite variation. 

5. MODELS WITH LEVEL DEPENDENT VOLATILITY 

We will display the method in the case of a complete market, the probability being 
a martingale measure and for a simple one-dimensional diffusion model. 
The stock is supposed to be the solution of the equation 

dX t = X t <j(t, X t ) dB t + X t r(t) dt. 

We limit the study to the error due to (Bt) which is defined by an Ornstein-Uhlen- 
beck structure: 

r[f °°h(s) dB s } = f °°h 2 (s)ds. 

The rate is deterministic, the function a(t,x) will be supposed bounded with 
bounded derivative in x uniformly for t G [0, T\. 

Let f(X T ) be a European option. Its value at time t is 

V t = neM-ft r r(s)ds)f(X T )\F t ] 

the hedging portfolio is given by the adapted process H t which satisfies 

(5.1) V t = exp(- S*r(s)ds)Vt = V + f* H s dX s 

where X s = exp(— J* * r(s)ds)X t . 

We proceed as follows: from the equation 



X t = X + f X s a(s,X s )dB s + fr(s)X, 
Jo Jo 



s ds 



we obtain 



X*= [ U (a(s,X s )+X s a' x (s,X s ))XfdB s + f ' X s a(s,X s ) dB s + f ' r(s)X#di 
Jo Jo Jo 

this equation is solved by putting 

K s = a(s,X s )+X s a' x (s,X s ) 

M u = exp {/ U K s dB s - \ / U Kids + J U r (s)ds} 



and remarking that 



Jo 



a) Let us first suppose / G C 1 fl Lip and let us define Y = exp(— J t T r(s)ds)f(X T ). 
To compute (!E[V|jF t ]) # we apply the second lemma of section 3: 

Y** =exp(-j t T r(s)ds)f (X T )Xp 



and 

iny\H)* = zm- si r{s)ds)nr{x T )xp\F t } 

= exp(- jf r{s)ds)nf{X T )M T \T t ] f X ° a( f > X ^ d B s 

Jo M s 

and the second lemma gives 

v[v t }= nnrim 

[ ] = exp(-2/ t T r( g )^)(E[r(X T )M T |^]) 2 /o ds 

this yields also the cross error of V s and Vt which is usefull to compute errors on 
random variables such that J T h(s)dV s or J T h(s)V s ds. 

T[V s ,V t ] = exp(- Jj 'r(s)ds - g r(s)ds) 
(5.3) nf{X T )M T \T s W{X T )M T \T t ] 

rsAt XZa 2 (u,X u ) i 
JO Ml 

With our hypotheses as t | T 

m) - nx T )M* £ x>2 ^ Xs) ds = nx T )v\x T ] 

in -^ 1 (1P) and a.s. 

b) Now to deal with H t , let us remark first that H t is easily obtained by the Clark 
formula. The formula 15 . II gives 

H t exp(-J Q t r(s)ds)X t a(X t ) = D ad [exp(- J T r(s)ds)f(X T )} 

where D ad is the adapted O-U-gradient defined by 

D ad [Z]{t) = ^[DZ{t)\T t ]. 

Since 

D[exp(-J T r(s)ds)f(X T )} = exp(- £ ' r(s)ds)f(X T )(DX T )(t) 
we have from the computation done for Vt 

D[exp(-gr(s)ds)f{X T )} = exp(- / T r{s)ds)M\f'(X T )M T \^ t ] 



Thus 



H t = exp(- J t T r(s)d S )nf'(XT)M T \F t }± 

Now supposing / and /' G C 1 fl Lip we apply the same method as for obtaining 
T[V t ] which leads to 

r[H t ] = exp(-2 /, r r(s)ds) 

(5.4) [ n ML (r[XT )M T + nx T) zrm) 2 f %? Hn - AV 



Ml 



du 



with Zj = / L s gLB s - / K s L s M s ds 
Jt Jt 

and K s = a(s,X s ) + X s a' x (s,X s ) 

L s = 2a' x {s,X s )+X s a^{s,X s ). 



If we introduce the following notation which, in our present Markovian model, 
gives the probabilistic interpretation of the usual Greeks 

delta, = H t = exp(- r(s)ds)1Eif' (X T )M T \F t ]± 
gamma, = exp(- f t T r(s)ds)n$(f"(X T ) + ^f{X T )Zj\F t ] 

we can summarize some formulas of this case with level dependent volatility by 



v* = 

m\ = 


delt&tXf 
delta t 2 rLY t ] 


m,v t ] = 


delta s delta t r[X s , X t ] 


rm = 


gamma t X t # 
gamma^r[X t ] 


T[H a ,H t ] = 
T[V s ,H t ] = 


gamma s gamma t r[X s , X t 
delta s gamma t r[X s , X t ] 


nx t ] = 


u 


r[x s ,x t ] = 


M s MtC tXW ^ Xu) du. 

1V1 U 



6. CONCLUDING REMARKS 



The error calculus based on Dirichlet forms begins at present to be used by mod- 
elisators in economics and finance. It is too early to give an account of its ap- 
plications. What we have done is just presentating how this tool runs through 
stochastic models including SDE's. 

Among the directions of research let us mention that this approach yields new 
integration by parts formulas which have been shown to be useful to compute the 
Greeks by Monte Carlo methods. Also, it allows to perform Malliavin calculus on 
the Monte Carlo sample space, that is after discretization instead of before. Let 
us indicate briefly the idea : 

Consider the error structure 

(n, A, P, ID, r) = ((0, 1), B(0, 1), dx, ctt, 7 ) N * 
dl= {u e L 2 (0, 1) x(l - x)u'(x) e H%(0, 1)} 
7 [«](#) = x 2 (l — x) 2 u' 2 (x). 

Let us denote U n the coordinate maps, taking H = £ 2 this structure admits the 
following gradient : if F — f(U 1: U 2 , ■ ■ ■ , U n , . . .) belongs to ID 

DF = (ti(U u U 2 ,...,U n ,.. )U n {\ - C/ n )) n >! 
and if a G £ 2 we have the following integration by parts formula 

E[< DF, a >fi] = -E[F£^ =0 a n (l - 2U n )]. 
Let us take for example the following discrete approximation of an SDE: 

S n +i = S n + a(S n )(Y n+ i — Y n ) , So = x, 

where Y n+ \ — Y n = \£(n + 1, U n+ i), we get easily under regularity assumptions on 
£ and a : 

£E[»(S W )]= wmwr u %&3g**-&fi 



As a second direction of research, let us sketch how to do a sensitivity analysis 
of the solution of an SDE with respect to a functional coefficient. 

Let us consider that the level dependent volatility of section 5 is a function in 
a vector space Q± equipped with an error structure (Qi, A\, Pi, IDi, I\) such that 
IPx-almost every function in Qi be of class C 1 and Lipschitz, and that the linear 
form V x defined by 

V x (a) = a(x) 

belong to Di. Defining the notation cr* by <r#(x) = (V x )#(a) we can show that if 
Y is an erroneous independent random variable defined on an other error structure, 
the following formula holds 

(a(Y))# = a*(Y) + a'(Y)Y*. 

Then similar computations to those of section 5 can be done. For example suppose 
a is represented for numerical evaluation as 

a(t,x) = ^2a pq t p x Q 
p,i 

where the CLpq S £1X6 random and erroneous with r[a P9 ] = aj^, T[a pq ,aij 1 = if 
(V) l) we obtain that the error on a transfers to X t in the following way: 

/ rt QPY q+l \ 2 

r[X t] = ^[j o ^-( d B,-K, is )) 4„ 

where M t , K t have the same meaning as in section 5. 

Because estimates of biases are important in financial models especially for 
pricing, see e.g. Hull and White (1988), let us mention shortly what would be 
the second order calculus with variances and biases mentionned above in the table 
2.1 of section 2. In an error structure (Q, A, IP, ID, T) the bias of the error on a 
random variable X (i.e. the conditional expectation of the error) is represented by 
the generator A of the semi-group canonically associated with the error structure 
acting on X, see Bouleau and Hirsch (1991). It has a domain DA smaller than ID. 
The functional calculus on A follows the following rules: for all F G C 2 (H d ), V/j 
locally in DA, i = 1, . . . , d, then F(fi, . . . , fa) is locally in DA and 

(6-1) A[F(f)} = y tn{f)Af i + \ E /•'/'(/)!' J,\ 

i=i z i,j=i 

On the Black-Scholes model, with the O-U hypotheses and concerning solely the 
error due to (B t ), we obtain : 

T[B t ] = t, T[S t ] = S 2 t a\ A[B t ] = -B t , 

A[S t ] = -StaB t + \a 2 S t t, A[V t ] = delta^] + |gamma t r[^], 

A[H t ] = gamma t A[S t ] + |g(t, S t , a, r)T[S t ]. 

Except at time t — (since the perturbation 3.3 doesn't move at t — but we 
can imagine it starts farther in the past) we see that an error on the path of the 
stock (5' t ) t >o induces biases on the price V t and on the hedge H t (involving the 
'Greek' §Jr)- May these biases due to a fuzzy stock price be an interpretation of the 



bid-ask ? This interesting question needs certainly more complete investigations 
since when a transaction occurs the price is erroneous while the amount of stock 
which is bought or sold is not. Anyhow it is an auspicious project to understand 
and to take in account the consequences of the bid-ask on pricing and hedging 
procedures of tractable financial models with the help of the error calculus on 
variances (operator T) and biases (operator A). 

REFERENCES 

Beurling, A., Deny, J. (1958-59): Espaces de Dirichlet, I. le cas element aire, 
Acta Math. 99, 203-224 ; Dirichlet spaces, Proc. Nat. Acad. Sci. U.S.A. 45, 
206-215. 

Bogachev, V. I., Roeckner, M. (1995): Mehler formula and capacities for in- 
finite dimensional Ornstein-Uhlenbeck processes with general linear drift, Osaka 
J. Math. 32, 237-274. 

Bouleau, N. (1995): Construction of Dirichlet structures, in : Potential theory 
ICPT 1994, de Gruyter. 

Bouleau, N. (2001): Calcul d'erreur complet lipschitzien et formes de Dirichlet, 
J. Math, pures et appl. 80, 9, 961-976. 

Bouleau, N., and Hirsch, F. (1986): Proprietes d'absolue continuity dans les 
espaces de Dirichlet et application aux EDS, in Sem. probabilite XX, Lecture 
notes in Mathematics 1204, Springer, 131-161. 

BOULEAU, N., and HlRSCH, F. (1991): Dirichlet forms and analysis on Wiener 
space, De Gruyter. 

Bouleau, N., and Lamberton, D. (1989): Residual risks and hedging strategies 
in Markovian markets, Stochast. Process. Appl. 33, 131-150. 

Dellacherie, and Cl., Meyer, P. A. (1987): Probabilites et potentiel. Her- 
mann. 

Feyel, D., and LA Pradelle, A. de (1989): Espaces de Sobolev Gaussiens, 

Ann. Inst. Fourier, 39-4, 875-908. 
Fournie, E., Lasry, J.-M., Lebuchoux, J., Lions, P.-L., and Touzi, N. 

(1999): Application of Malliavin calculus to Monte Carlo methods in finance, 

Finance and Stochastics, 391-412. 
FUKUSHIMA, M., OSHIMA, Y., and TAKEDA, M. (1994): Dirichlet forms and 

Markov processes. De Gruyter. 
Hull, J. and White, A. (1988) An analysis of the bias in option pricing caused 

by a stochastic volatility, Advances in Futures and Options Research, vol 3, 

29-61. 

Lamberton, D., and Lapeyre, B. (1997): Introduction au calcul stochastic 

applique a la finance. Ellipses. 
Ma, Z., and ROECKNER, M. ( 1992): Dirichlet forms. Springer. 
MALLIAVIN, P. (1997): Stochastic analysis. Springer. 

NUALART, N. (1995): The Malliavin calculus and related topics. Springer. 



